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Abstract 

We obtained accurate resonances for the Stark effect in hydrogen by means of three 
independent methods. Two of them are based on complex rotation of the coordinates 
and diagonalization of the Hamiltonian matrix (CRLM and CRCH). The other one 
is based on the Riccati equations for the logarithmic derivatives of factors of the 
wavefunction (RPM). The latter approach enabled us to obtain the most accurate 
results and extremely sharp resonances. 


1 Introduction 


The Stark effect in hydrogen is an old problem in atomic spectroscopy and 
one of the first triumphs of wave mechanics [HE] (and references therein). 
The Schrodinger equation is separable in parabolic and squared parabolic 
coordinates which facilitates the application of most approximate methods [2|- 

In a recent paper Fernandez-Menchero and Summers [3] obtained the complex 
eigenvalues and eigenfunctions of the Hamiltonian operator for the hydrogen 
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atom in a uniform electric field. They resorted to the Lagrange-mesh basis set, 
proposed by Lin and Ho [1] for the treatment of the Yukawa potential in a 
uniform electric held, and the complex-rotation method [S]. They compared 
their results with those obtained by Lin and Ho |1], Kolosov |6], Rao and Li 
and Ivanov [8] and overlooked the earlier impressive calculations of Benassi 
and Grecchi [S] and the accurate results obtained by Fernandez PI- Benassi 
and Grecchi resorted to complex scaling and a basis set of conhuent hyper¬ 
geometric functions that is suitable when the Schrodinger equation is written 
in squared parabolic coordinates. On the other hand, Fernandez applied the 
straightforward Riccati-Pade method (RPM) that does not require the use of 
complex coordinates. 

The purpose of this paper is to calculate the Stark resonances as accurately 
as possible by means of the methods proposed by Fernandez-Menchero and 
Summers [3], Benassi and Grecchi [9] and Fernandez [10] and compare the 
results with those obtained by the authors already mentioned and also by 
Hamburg and Kolosov HU. There is a vast literature on the hydrogen atom 
in a uniform electric held but we restrict present discussion to some of the 
available calculations that we deem are suitable for comparison. 

In section [2] we outline the main ideas about separating the Schrodinger equa¬ 
tion in parabolic and squared parabolic coordinates. In sections [31 ID and [3] 
we briehy introduce the methods of Fernandez-Menchero and Summers [3], 
Benassi and Grecchi [9] and the RPM [10], respectively. In section [6] we com¬ 
pare the results of various approaches and in section [7| we summarize the main 
results and draw conclusions. 


2 Stark effect in hydrogen 


The Schrodinger equation in atomic units is 
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( 1 ) 


H=-l-V^-- + Fz, 

2 r 

where F is the intensity of the uniform electric held assumed to be directed 
along the 2 ; axis. 

This equation is separable in parabolic coordinates 



^ > 0, > 0, 0 < 0 < 27r. (2) 

If we write 

y, z) = m = 0, ±1, ±2,..., (3) 


then we obtain two equations of the form 


f cP 1-rrP E_ F 

\dx‘^ 4a;^ 2 A ^ x j 


<I>(a:) = 0, 


(4) 


where a = ±1 and ^ 4 + = A and = 1 — ^4 are separation constants. When 
(T = 1, X = (^ and when a = —1, x = rj and $( 77 ) = v{ri). 

The Schrodinger equation ([T]) is also separable in squared parabolic coordi¬ 
nates 


X = y.i'coscj), y = y.i'smcj), 



/i>0, i' > 0, 0<(/)< 27r. 


If in this case we write 
P{x,y,z) = 


(5) 


( 6 ) 
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then we obtain two equations of the form 


+ 2Ex^ - aFx^ + Z,j <h(x) = 0, (7) 

where, a = ±1 and Z_|_ = Z and Z_ = 4 — Z are the separation constants. 
When cr = 1, a: = yU and d’(/i) = 'u(/i); when a = —1, x = u and $( 1 ^) = v^u). 

The solutions to the equations in either set of coordinates are commonly la¬ 
belled by the quantum numbers ni,n 2 = 0, 1 , 2 ... and m = 0, 1 ,.. and the 
notation \ni,n 2 ,rn) is suitable for referring to them. We will sometimes resort 
to the principal quantum number n = ni -|- n 2 -|- |m| -|- 1 to denote a set of 
states. Obviously, m is the only good quantum number; the other ones refer 
to the states of the hydrogen atom and are valid when F = 0. 


I (P 1 — 4m^ 
\dx‘^ 4a;2 


3 Complex rotation and Lagnerre-mesh basis set 


Fernandez-Menchero and Summers [3] decided to treat the Schrodinger equa¬ 
tion as nonseparable. The Hamiltonian operator in parabolic coordinates reads 

\ d f d 

j + 

and the authors proposed the variational ansatz 
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2 ^rj dcfP i + Tj 


+ F 


^-v 


,( 8 ) 


H = 


^ + V 


di 'Cdi 


N N 


CklmG ^ 2 ’' ANkiOPmiv) 

y Zn 


ANk{x) = 


k=l 1=1 

Ln{x) 


X-Xk 


(9) 

( 10 ) 


where L]^{x) is the Laguerre polynomial of degree N and Xk its A;-th zero. In or¬ 
der to obtain the resonances they resorted to the well-known complex rotation 
method [5] that in this case is given by the transformation (^, rj) —)■ 
where "d is the rotation angle. The eigenvalues and expansion coefficients are 
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given by the secular equation 


{U-ES)C = 0, 


( 11 ) 


where the elements of the x N‘^ matrices H and S are explicitly shown 
elsewhere |3] and the elements of the column vector C are the coefficients 


Ckim- Note that the integrals appearing in the matrix elements of both H 
and S should be calculated numerically and when we increase N we have to 
calculate all those integrals again. For brevity we will call this method CRLM. 


4 Complex scaling and confluent hypergeometric basis set 


In order to obtain the resonances Benassi and Grecchi [5] resorted to equation 
([7]) and a basis set of the form 



( 12 ) 


where F{a, b, z) is the confluent hypergeometric function. In this case the 
authors resorted to the complex scaling method that is based on the transfor¬ 


mation —)■ where A is a complex number. The complex 


scaling method contains the complex rotation method as a particular case 
because A = |A|e*’^ and the proper choice of |A| enables one to improve the 
convergence of the approach. 

In this case all the elements of the relevant pentadiagonal matrix can be 
calculated analytically and are independent of the matrix dimension. This 
approach will be called CRCH from now on. 
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5 The Riccati-Pade method 


We can apply the RPM to the eigenvalue equations derived in either parabolic 
or squared parabolic coordinates. In the earlier application of the approach 
Fernandez HDl chose the former and here we resort to the latter. It is worth 
mentioning that the performance of the RPM in both sets of coordinates 
is identical and that the reason for the selection of the squared parabolic 
coordinates is to have a closer contact between the RPM and the CRCH 
method of Benassi and Grecchi [5]. The regularized logarithmic derivative 


/(^) 


s ^'(x) 
X <h(x) ’ 


1 

2 ’ 


(13) 


can be expanded in a Taylor series 

OO 

fix) = (14) 

j=o 


where the coefficients fj are polynomial functions of E and Z. The details of 
the method are outlined elsewhere [lO]; here it suffices to say that we construct 
Hankel determinants of the form 


hUe,z,f) = 


fd-\-\ fd-\-2 


fn+d 


fd+2 fd+3 ■ ■ ■ fn+d+l 


(15) 


fn+d fn+d+l ■ ■ ■ f2D+d-l 


and obtain the approximate eigenvalues from the roots of the set of 

nonlinear equations 

Hi{E, Z, E) = HUE, 4 - Z, -E) = 0. (16) 
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The main advantage of the RPM is the enormous rate of convergence which 
enables us to obtain very accurate eigenvalues with determinants of relatively 
small dimension. However, the great number of roots in the neighborhood of 
each eigenvalue makes it difficult to hnd the optimal sequence that converges 
to it. Since we resort to the Newton-Raphson algorithm to obtain the roots of 
the system of equations flTBD we have to choose the starting point quite close 
to the chosen root. We will discuss this point briefly in Section [HI The RPM 
is most suitable for the treatment of separable problems. 


6 Numerical calculations 


In order to apply the CRLM [3] we calculated all the integrals numerically with 
a tolerance of 10“^®. For each value of F, m and N we varied the rotation angle 
■d between 0.3 and 0.7 looking for those eigenvalues that remained almost 
constant. We could reproduce the results in the literature with matrices of 
dimension = 30 [12] • 

The only exact quantum number is m, however, it is customary to resort to 
the quantum numbers of the isolated hydrogen atom in order to label the 
energies and states of the Stark problem. Some authors choose the parabolic 
quantum numbers ni,n 2 = 0,1,2,... |9] and others the principal quantum 
number n = rii + n 2 + |m,| -|- 1 and k = rii — n 2 [5]. 

In the case of CRCH we hrst solved the eigenvalue equations that yield the 
eigenvalues and Z_ and then applied the Newton-Raphson method to 
solve the equation -I- Z_ — 4 = 0. Details of the calculation are given 
elsewhere [9]; we just mention that in order to obtain a starting point for the 
Newton-Raphson method we resorted to the results provided by perturbation 
theory |l6j . 
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The Hankel determinants that appear in the RPM are polynomial functions 
of E and Z of great degree. For this reason it is necessary to handle complex 
numbers with sufficiently great precision and we resorted to the GNU MPC 
library [13] . The Hankel determinants can be calculated numerically by means 
of the well known recurrence relation 



with the initial conditions Hq = 1 and Hf = fd+i- Once we calculate the 
desired determinants we obtain the eigenvalue and separation constant by 
means of the Newton-Raphson method. In order to have a suitable starting 
point we resorted to CRCH results. 


The remarkable rate of convergence of the RPM is clearly illustrated by the 


calculation of the logarithmic error log 


where a stands for either 


alD] _ 

ReU or \m.E. We do not indicate the value of d explicitly because in present 
calculations we have chosen d = 0. Figures [D and [2] show the logarithmic 
error for all the resonances with n = 1,2,... ,6 and field strengths F = 0.001 
and F = 0.005, respectively. As the quantum numbers increase the minimum 
value of D at which the resonance appears also increases. For example, the 
lowest resonance appears as a root of the Hankel determinant with D = 2 and 
|0, 5, 0) appears at D = 11. The rate of convergence is greater when the root of 
the Hankel determinant is real. This fact is clearly shown in Fig. H] where the 
imaginary part of |0, 0, 0) appears at D = 22. When F = 0.001 the imaginary 
part of the lowest resonance appears at D = 103 which is the reason why the 
rate of convergence for this resonance is considerably greater than the other 
ones for all D < 100. When D > 103 the rate of convergence for the lowest 
resonance becomes similar to the other ones. 


The resonances explicitly labelled in figures [T] and [2] are shown in Tabled] with 
their number of digits truncated to a reasonable size. 





Figure |3] shows that our estimated value of \m.E for the lowest resonance is 
in perfect agreement with the analytic asymptotic formula derived by Benassi 
and Grecchi [S]: 


|ImE| ~ _ 8.916F + 25.57^2 + 0{F^)) . (18) 


The RPM is also suitable for the calculation of higher resonances. For example, 
tables [2] and [3] compare present results obtained by the CRCH and RPM for 
two states with n = 5 with those obtained earlier by Damburg and Kolosov 
[m. Table m compares present RPM results for some states with n = 10 with 
those obtained by Kolosov j6]. The discrepancy in the imaginary part for the 
case |0, 9, 0), F = 2.2 x 10“^, is probably due to a misprint in that reference. 

Table 0 shows the resonance |39,0,0) for several values of the held strength. 

We do not compare these results with those of Kolosov [6] because he did not 
indicate the conversion factor from atomic units to Vcm~^ shown in his table. 

However, Fig. 0] shows that both sets of results are in reasonable agreement. 

In the tables discussed above we have truncated present RPM results to a rea¬ 
sonable number of digits. We have obtained them with much higher accuracy 
as suggested by hgures [T] and [2l For example, for the lowest resonance and 
held strength F = 0.005 we obtained 

Re^ =-0.50005628479379296933177394769143288196325092731889137262135731 
28725736315548994436307340293823812601699152241599625041068943791 
42099665225189334039046848974164185728077545219665133771938893895 
64251327341968732189236225621425838831553440690618168917215735013 
803880912033036 

Im^=-4.74901370837102040886757127120827250441845432417751748825418970 
22400488040285011762035775189238632536585799373474503067879411046 
22147574080708907330396144467615023762954201754322979890803189455 
51562966634796276868224 x 10“^® (19) 

with Hankel determinants of dimension D < 150. In principle we expect that a 
properly truncated perturbation series will exhibit an accuracy of the order of 
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|Im£'|. On summing the first 130 terms of the perturbation series calculated by 
means of the hypervirial perturbative method [16] we obtained the following 
result: 


= -0.5000562847937929693317739476914328819632509273188913726, (20) 


that agrees with the RPM one to the last digit. It is not easy to obtain such 
a sharp resonance by means of other approaches; for example, Fernandez- 
Menchero and Summers [3] estimated Rei7 = —0.5000553416 and Im£^ = 
0.8944475605 x 10“^. We calculated the real part of this resonance more accu¬ 
rately by means of the CRLM but were unable to obtain a reasonable estimate 
of the imaginary part |12] . 

The RPM enables one to calculate even sharper resonances; for example, for 
the lowest one and F = 0.001 we obtained 


ReF;= -0.50000225005555178356591589970608204532866714376652965654995937 
97019283545891048870035463753481536961447568150634794700138591827 
91549628581187487453336046428670620173909589867079695807271725700 
47474205292728633151353049600188535220623998127315129221076077663 
756392409425470889188167975544640438386213612059475282765271923 
ImE = -5.854592875137598393486482622915575 x 10”^®^ (21) 


with Hankel determinants of dimension D < 130. In this case perturbation 
theory of order 600 (300 nonzero terms) yields 


=-0.50000225005555178356591589970608204532866714376652965654995937970 
19283545891048870035463753481536961447568150634794700138591827915496 
28581187487453336046428670620173909589867079695807271725700474742052 
92728633151353049600188535220623998127315129221076077663756392409425 
47088918816797549 (22) 
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7 Conclusions 


We have calculated the resonances of the Stark effect in hydrogen by means 
of three independent methods. Although we were able to improve the CRLM 
calculation considerably |T2] we think that the CRCH is far more efficient. 
However, the RPM yielded considerable more accurate results and enabled 
us to obtain extremely sharp resonances that we were not able to obtain by 
means of the other two methods. The reason is that the accuracy of the real 
part should be at least of the order of the imaginary one. We were able to attain 
such an accuracy in the calculation of the roots of the RPM equations fll6l) 
thanks to the GNU MPC library [13]. We think that it is almost impossible 
to do the same by means of the CRLM because of the numerical calculation 
of the matrix elements. In principle, one can obtain the resonances with any 
degree of accuracy by means of the CRCH but such calculation would require a 
great deal of ingenuity. For this reason we think that the RPM is an extremely 
suitable benchmark to test other approaches on separable models. 
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Table 1 


Resonances for the states appearing in figures [U and [2] 


F 

Resonance 

Re{E) 

Im(F;) 


|0,0,0) 

-0.5000022500555518 

-6.584169959231863 x 10"287 


|0,0,1) 

-0.1250782240371032 

-8.433615180808857 x lO'^s 

0.001 

|0,1,0) 

-0.1280858350607099 

-2.060525710039887 x lO'^i 


|1,0,0) 

-0.1220826861326878 

-3.395926205766083 x lO'^^ 


|0,5,0) 

-0.05215538955477732 

-2.594493723108199 x lO"^ 


|0,0,0) 

-0.5000562847937930 

-4.749013708371020 x 10"®® 


|0,0,1) 

-0.1271466127039709 

-1.307642723230557 x 10"^ 

0.005 

|0,1,0) 

-0.1426186075727077 

-5.297223183652474 x 10"^ 


|1,0,0) 

-0.1120619240019938 

-2.864684219868783 x 10"® 


|0,5,0) 

-0.1213596730003857 

-1.176260968442979 x 10"^ 


Fig. 1. Convergence of the RPM resonances for F = 0.001 
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Table 2 

Resonance |0,4,0) from reference [TT] (a) and present calculation by means of CRCH 


(b) and RPM (c) 


F 


Re{E) 

ImlE) 


a 

-0.0231791962 

-2.1 X 10-12 

0.00010 

b 

-0.02317919625030 

-2.1135 X 10-12 


c 

-0.02317919625030518 

-2.113884073268850 x 10-12 


a 

-0.024956749 

-9.595 X 10-1 

0.00015 

b 

-0.024956750918078 

-9.6007202913 x IQ-i 


c 

-0.02495675091807878 

-9.600720291331372 x IQ-i 


a 

-0.02697136 

-8.9150 X 10-® 

0.00020 

b 

-0.0269800814710915 

-9.36280360832384 x IQ-^ 


c 

-0.02698008147109154 

-9.362803608323849 x IQ-^ 


a 

-0.02896828 

-4.2655 X 10-1 

0.00025 

b 

-0.02912946983310681 

-4.868994650393436 x IQ-i 


c 

-0.02912946983310681 

-4.868994650393436 x IQ-i 


a 

-0.0305381 

-9.849 X 10-1 

0.00030 

b 

-0.03122955458572655 

-1.127494087615666 x IQ-^ 


c 

-0.03122955458572655 

-1.127494087615666 x IQ-^ 


a 

-0.0314338 

-1.8217 X 10-3 

0.00035 

b 

-0.03323652729915596 

-1.945324601526496 x IQ-^ 


c 

-0.03323652729915596 

-1.945324601526496 x IQ-^ 


a 

-0.031408 

-3.17565 X 10-3 

0.00040 

b 

-0.03512209724011620 

-2.892832253491630 x 10-3 


c 

-0.03512209724011620 

-2.892832253491630 x 10-3 


a 

-0.02998 

-6.365 X 10-3 

0.00045 

b 

-0.03687445248862566 

-3.911655467856354 x 10-3 


c 

-0.03687445248862566 

-3.911655467856354 x 10-3 
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Table 3 

Resonances for the state |4,0,0) from reference [11] (a) and present calculation by 


means of CRCH (b) (b) and RPM (c) 


F 


Re{E) 

ImlE) 


a 

-0.0158077645 

-1 X 10-11 

0.00015 

b 

-0.01580776440749585 

-7.156147028941416 x lO^i^ 


c 

-0.01580776440749585 

-7.156147028941416 x lO^i^ 


a 

-0.0145352049 

-2.013-8 

0.00020 

b 

-0.01453520517676726 

-2.012419057345574 x lO^® 


c 

-0.01453520517676726 

-2.012419057345574 x lO^® 


a 

-0.013328925 

-1.63595 X 10-® 

0.00025 

b 

-0.01332892813256598 

-1.637235677233378 x lO^® 


c 

-0.01332892813256598 

-1.637235677233378 x lO^® 


a 

-0.01220093 

-2.0833 X 10-® 

0.00030 

b 

-0.01220135935615766 

-2.104916128829678 x lO^® 


c 

-0.01220135935615766 

-2.104916128829678 x lO^® 


a 

-0.01113604 

-8.9570 X 10-® 

0.00035 

b 

-0.01114288854595917 

-9.327043407081445 x 10^® 


c 

-0.01114288854595917 

-9.327043407081445 x 10^® 


a 

-0.01008206 

-2.1402 X 10-1 

0.00040 

b 

-0.01011729953739499 

-2.321591626792999 x lO^i 


c 

-0.01011729953739499 

-2.321591626792999 x lO^i 


a 

-0.00899479 

-3.7941 X 10-1 

0.00045 

b 

-0.00909725070184054 

-4.263615594700631 x IQ-i 


c 

-0.00909725070184054 

-4.263615594700631 x l-i 


a 

-0.0078517 

-5.7415 X 10-1 

0.00050 

b 

-0.00807076238659657 

-6.601708601710509 x IQ-i 


c 

-0.00807076238659657 

-6.601708601710509 x IQ-i 
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Table 4 

Resonances calculated by Kolosov [6] (a) and by means of the RPM with D = 30 

(b) 

|9,0,0) 


10® F 


ReF 

ImF 


a 

-2.58557398 x 10^® 

-9.509227 X 10^* 

2.0 

b 

-2.585573979364734 x lO^® 

-9.509226978683682 x lO^® 


a 

-1.57105982 X 10^® 

-3.959433 x lO^® 

3.0 

b 

-1.571059822031523 x 10^® 

-3.959432995212450 x 10^® 


a 

-5.8496223 x 10^'^ 

-1.6703408 X 10-^ 

4.0 

b 

-5.849621042229387 x lO^'^ 

-1.670340346132577 x 10^^ 


|0,0,9) 


10® F 


ReF 

ImF 


a 

-5.32440479 x 10^® 

-3.9351431 X 10-® 

2.0 

b 

-5.324404794258087 x 10^® 

-3.935143048784509 x 10^® 


a 

-5.6483507 x lO^® 

-3.263613 X 10-^ 

3.0 

b 

-5.648350339949772 x 10^® 

-3.263623549768137 x lO^'^ 


|0,9,0) 


10® F 


ReF 

ImF 


a 

-7.2120845 X 10^® 

-4.0070215 X 10-® 

1.4 

b 

-7.212084472616482 x 10^® 

-4.007021552503371 x 10^® 


a 

-7.977367 x lO^® 

-2.4391785 x lO^'^ 

1.8 

b 

-7.977367228278029 x 10^® 

-2.439179991742931 x 10^^ 


a 

-8.660578 X lO^® 

-5.32992 X 10-® 

2.2 

b 

-8.660579416493959 x 10^® 

-5.329919686471733 x 10^^ 
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Table 5 


Resonance |39,0,0) calculated by means of the RPM with D < 65 


lO’^F 

ReS 

Im E 

1.0 

-1.033131815036742 x 10“^ 

-1.410563208376528 x lO'^^ 

1.2 

-6.597779434524293 x lO"'^ 

-2.882918724695370 x 10"® 

1.4 

-3.007626411426787 x 10"^ 

-8.858397640808244 x 10"^ 

1.6 

5.730328096956075 x 10"® 

-2.598217010741238 x 10"® 

1.8 

4.153413722356150 x 10“^ 

-4.277370190378267 x 10"® 

2.0 

7.701721282910781 x 10"^ 

-5.720318291587286 x 10"® 

2.2 

1.120415206875958 x lO'^^ 

-6.975702947269457 x 10"® 


Fig. 2. Convergence of the RPM resonances for F = 0.005 





Fig. 3. Width of the lowest resonance calculated by means of the RPM (circles) and 
the asymptotic expansion (fTSll (line) 
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Fig. 4. Real and imaginary parts of the resonance |39,0,0): (a) [6], (b) RPM 
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